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I. INTRODUCTION 



In a recent paper a multichannel algebraic scattering (MCAS) theory for nucleons 
scattering from a nucleus was specified in detail. At low energies, the approach is noteworthy 
because its formulation facilitates a systematic determination both of the sub-threshold 
bound states and of the compound resonances. The theory is built upon sturmian expansions 
of whatever one chooses to be an interaction matrix of potential functions Of course there 
is the usual limitation of dealing with a coupled channels problem namely that the dimension 
of the evaluation rapidly increases with the number of channels considered. However, the 
MCAS approach does treat all selected channels equivalently, whether they be open or closed, 
so that solutions can be found at both positive (scattering) and negative (bound system) 
energies relative to the nucleon on nucleus threshold. 

While the MCAS approach may be used for any target (and projectile) system, our 
formulation to date has been for nucleon-nucleus interactions. However, the most practical 
adaptations insofar as size of the problem and, concomitantly, computation times, are with 
light mass targets. They have well separated low excitation spectra and usually the total 
cross-section data show distinct resonances upon a smooth background. Thus the first 
MCAS study P was of low-energy neutron- 12 C (n+ 12 C) scattering. For that only three 
states of 12 C ostensibly were needed in the evaluations; namely the 0^ (ground), l\ (4.4389 
MeV), and the (7.6542 MeV). Herein we consider that same system in more detail and, 
as well, analyze proton- 12 C (p+ 12 C) data. 

The MCAS method solves the coupled-channel Lippmann-Schwinger equations for the 
nucleon-nucleus system considered. The starting matrix of potentials may be constructed 
from any nuclear model: shell or collective, rotational or vibrational, model. We have used 
a rotational collective-model representation for the interaction potential matrix but therein 
take deformation to second order. We have chosen Woods-Saxon functions and their various 
derivatives to be the form factors for all components each with characteristic operators of 
central, spin-orbit (1-s), orbit-orbit (1-1), and of spin-spin (s-I) type. The interactions were 
allowed to depend on parity as well. With such a characterization, the model potential is 
sufficiently flexible to describe all possible structures (at both positive and negative energies) 
of the nucleon- 12 C system. However any such collective model prescription violates the Pauli 
principle. At these energies that violation is severe. But it is possible 0, 0| in the MCAS 
approach to account for Pauli blocking of occupied nucleon states in the target, with this or 
any collective model specification of the matrix of potentials. That is achieved by introducing 
orthogonalizing pseudo-potentials 0, 0, HI (OPP) into the scheme by which the sturmians 
are specified. In that way all sturmians in the (finite) set selected as the basis of expansion 
of the matrix of potentials contain few or no components equivalent to the external nucleon 
being placed in an already densely occupied orbit. Treatment of the Pauli principle has a 
significant effect on results []], |j| . 

The aim of this paper is to extend calculations from the 13 C to the 13 N system, and to 
carry out a comparative analysis of bound and resonance spectra involved. A second version 
of the program has been set up, able to deal on the same footing with the scattering of 
neutrons and protons from the target nucleus. The neutron results confirm those of Ref. 

It will be shown that, in the n+ 12 C process, the spectrum of resonances up to about 
6 MeV (in the laboratory (lab) system) is almost completely described by a mechanism 
involving excitation of the first 2 + level of 12 C with energy €2 = 4.4389 MeV. The spectrum 
shows a sequence of compound resonances, generated by the ~ , | + , and | + bound (single 



2 



TABLE I: n+ 12 C potential parameters (strengths in MeV). 



parity 


central 


orbit-orbit 


spin-orbit 


spin-spin 




-49.1437 


4.5588 


7.3836 


-4.7700 


+ 


-47.5627 


0.6098 


9.1760 


-0.0520 


other parameters: 












tq = 1.35 fm. 


ao = 0.65 fm. 


/3 2 = -0.52 





nucleon) states in 13 C. The situation is very similar in the p+ 12 C process, with one overall 
energy shift due to the Coulomb interaction. Because of this shift in energy, some of the 
compound resonances in the n+ 12 C system become quasi-compound ones in the p+ 12 C case. 

We will show that when deformation (which for the nucleon- 12 C system we consider is 
of quadrupole type and so linked to a parameter /3 2 ) tends to zero, the compound and 
quasi- compound states tend to pure states. For compound resonances, their widths tend to 
zero and their centroid energies tend to those of single particle bound states plus the core 
excitations, e^. For quasi- compound resonances, the widths tend to the natural widths of 
the single particle resonance upon which they are formed, while their centroid energies tend 
to those of the single particle resonance plus the core excitations 6j. It is interesting then to 
analyze the behavior of the phenomenology contained in the model as (3 2 varies continuously 
from the deemed physical value to zero. There is a double purpose to this, namely to check 
the rules outlined above in a significant physical case and to describe the spectroscopy of 
13 C and 13 N in terms of the scheme specified. 

In Sec. UU we discuss the results on the n+ 12 C system, some of which were reported 
previously 1[. But now we analyze the origin of each state by studying results in the zero 
deformation limit. In Sec. 11111 the same analysis is extended to a study of the p+ 12 C system, 
where the Coulomb shift transforms some of the compound resonances into quasi-compound 
ones. All calculations of this system have been carried out assuming charge symmetry. The 
potentials that gave the results discussed in Sec. HH (and previously |l|) have been used 
without change. We present conclusions we have drawn from these studies in Sec. IIVI 



II. THE 13 C SYSTEM 

Calculations of the n+ 12 C system have been carried out with the parameters defining the 
initiating nucleon- nucleus interaction matrix of potentials being those used previously [l|. 
They are presented again in Table |l] for easy reference. Couplings of the input channel with 
the ground state {e x = 0), with the 2+ (e 2 = 4.4389 MeV) and with the 0^ (e 3 = 7.6542 
MeV) excited states of the 12 C target have been considered. Moreover, the Pauli principle 
has been taken into account throughout using the OPP procedure However, note that 

the OPP treatment we use is an approximation. We have used a large but finite strength 
in place of the infinite value that the OPP theory presumes. One must thus expect some 
small spuriosity, which may affect our calculated energies but by about a tenth of an MeV. 

In Fig.Q our calculated n+ 12 C elastic cross section is compared with the evaluated nuclear 
data file (ENDF) formed by Pearlman [fj. Source data and references were obtained using 
the computer index of neutron data (CINDA) 0]. The energies are in the laboratory frame. 
The sequence of resonances evident therein as energy increases have spin-parity assignments 
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FIG. 1: (Color online) Comparison between data (filled circles) and theoretical (solid curve) cal- 
culated elastic scattering cross sections for n+ 12 C scattering to ~ 5 MeV. 



| + , | + , | , | + , | + , | + , | + . Details of those resonances and of bound states are compared 



1]. 



with experimental values ^Sj in Table ITT1 We have found the parameters of the resonances by 
studying the trajectories of the sturmian eigenvalues using the algorithm defined in Ref. 1| 
This we defined as the resonance identification (RI) process. As already noticed in Ref. 
this choice is unambiguous in the case of narrow resonances, while for wider resonances it 
represents just one of the possible methods for resonances identification. Note that in this 
table, the energies of the states (in MeV) are given in the center of mass (cm) system. The 
widths, on the other hand, are expressed in keV. In the first column of Table |H] we give an 
identifying number (n) to each state which we use throughout the following discussion. 

In the spectrum, the | (n = 2) resonance is not seen experimentally. However, from our 
model calculation it is practically coincident with the first and strong | + resonance (n = 6), 
and is shown in the figure like a narrow spike over the broader | + peak. Possibly that is why 
it has not seen in experiments to date. But an isobaric analogue to this has been observed in 
13 N, and its nature is well explained within the general logic of the spectroscopic structure 
we believe to be sensible. The very narrow | resonance (n = 8) which lies just above 
threshold has a partner in our calculations, but as a bound state just below threshold. 

The peak shown in Fig. ^ at the energy of about 4.40 MeV in the laboratory frame 
(4.06 MeV in the cm system) corresponds to our calculated | + (n = 11) resonance. In the 

neighborhood the model also predicts (by using the RI procedure) a | + resonance (n = 4) 
whose centroid energy is E\ ah = 5.05 MeV (E cm = 4.66 MeV). No corresponding peak in the 
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TABLE II: Sub-threshold bound states and low energy resonances of the n+ C system. 
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0.745 


elastic cross 


section for this 


is seen. '. 


In fact, that resonance 


is missing entirely 


in Fig. ^ To 



solve this puzzle we calculated the even-parity components of the cross section separately; i.e. 
the cross section was calculated for the unique J n value of | + , then with | + , etc. The results 
shown in Fig. El clearly identify the component responsible for each even-parity resonance. 
The only peak not accounted for in the analysis is that having a spin-parity | . In particular 
under the maximum at about 4.4 MeV, there are two overlapping resonances. They are the 
| + component which has a sharp maximum in the energy region under consideration and 

the | + component revealed as a very small bump at about the same energy. It is a very 

broad resonance. The | + component is more easily recognized in the calculated cross section 
when, as we discuss later, the deformation is decreased. Thus the cross section is, as could 
be expected, a dominantly s-wave background with strong narrow and broad resonances 
superimposed. The individual values of energy at which the s-wave cross section initially was 
calculated are represented by the x in the plot for J n = | + . The clustering of x's (around 
regions in which sharp resonances in the total cross section occur) reflects the density of 
mesh points used to find the variation through those regions as precisely as possible. Such 
requirement is further evidence of the need in any such study to ensure that all resonance 
centroids and half widths are defined in the process, and energy steps in that region are 
selected appropriately, otherwise with too large an energy step they will be missed. The RI 
procedure |l| ensures all resonances will be found in the overall energy regime to be studied. 
Also built into the code implementing this procedure is an automatic increase in numbers of 
energy points. The span of those points is dependent upon the half width of the resonance, 
as defined by the RI procedure. 

We conclude therefore that the 4.4 MeV peak in the measured cross section is mainly a 
| + resonance. The ^ + resonance foreseen by the model calculations is weak and completely 
masked by it. The effect may be understood better by studying the (cm) energy variation 
of the i + (total) scattering phase shift. That variation is displayed in Fig. |21 This phase 
shift does not exhibit a regular increment between 0° and 180° (as do the other resonant 
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FIG. 2: (Color online) Components in the n+ 12 C elastic cross section. The theoretical curve (solid 
curve) is decomposed in its even-parity components. The ^ + component is shown by the generally 
smooth curve connecting ' x ' marks. It constitutes virtually the entire cross section at energies near 

o -f- c -f- 

threshold. The ^ and | components are portrayed by the dashed and dotted curves respectively. 
The narrow | + and | + resonances at (lab) energies of 2.84 and 4.88 MeV respectively are portrayed 
by long dashed curves. 



phase shifts), but start from 180° (because of Levinson's theorem) at zero energy, reaching 
a minimum of 84.5° at 3.17 MeV and then increases passing through 173° at 5.063MeV. 
The result is a broad and not very pronounced peak in the partial cross section, as shown in 
Fig. |21 Hence we assign our | + (n = 11) level with the well defined experimental resonance 

at E cm = 3.9314 MeV rather than as the | + state heretofore assumed in the literature 0]. 

The agreement between theory and experiment shown in Fig. ^ and in Table |H] is very 
good. Energies of both bound and resonance states, and the widths of those resonances, 
are quite well matched. Particularly so when one recalls that the target structure has been 
taken within a simple rotational model scheme. A measure of the goodness of fit is the mean 
square error in the table entries. That measure formed using 11 of the states is 

H = ^\/^(£th-£cx P ) 2 = °- 0776 MeV - W 

Entries n = 2 and n = 4 have been omitted because the first has been not detected experi- 
mentally while for the second, we propose a different assignment of quantum numbers. 

Our procedure enables us to interpret the structure of the spectra by following the general 
discussion that has been given [^] about the properties of compound resonances. Consider 
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FIG. 3: Energy variation of the calculated ^ 



2 scattering phase shift. 



first the even-parity states. These can be specified in the chosen representation in terms 
of pure states that are identified by the set of quantum numbers, {J n ,I,j, £}. While the 
relevant details concerning these sets of quantum numbers have been extensively discussed 
before 0, it is useful to note here that the single particle quantum numbers (/) and (J) can 
be related to a state of the single nucleon in the mass-13 systems. Then, in the energy range 
we study, and considering only the effect of the most important couplings, i.e. with the 
ground and 2 + state of 12 C, all possible configurations are as listed in Table ITTT1 The second 
excited 0^ state will enter discussion later in relation to a particular (| ) excited state in 
the mass-13 spectra. Assuming a shell model single nucleon spectrum for an additional 
nucleon coupling to states in 12 C we expect at most three even-parity bound states to be 
important, namely relating to lodgment of that nucleon within the \s-0d shell. However, as 
reported shell model studies suggest that the low excitation mass-13 states dominantly 
are identifiable with either Ods or lsi orbit couplings to the ground and 2 + states in 12 C. 

2 2 

As such, the coupling to the ground state (in 12 C) will provide two states in the compound 
nucleus with spin-parity values of | + and | + . Therefore we expect one doublet and one 
quintuplet of even-parity states when the coupling is made with the 2 + state. 

Using this scheme, the bound state | + (n = 3 of Table |TTJ) in 13 C should be dominantly 
described as a lsi neutron bound to the 12 C (0 + ) ground state core and identified with the 
N = 1 entry in the table of pure states: this infers a bound state of energy E ~ —2 MeV. 
Likewise there should be a | + state (n=9 of Table ITT]) . which can be identified with the N 
= 9 entry in Table IIHI Then, when a neutron impinging with energy 62 + E loses 62 to 
excitation and is bound to the 12 C* (2 + ) core, there should arise a very narrow resonance; 
a bound state in the continuum for zero coupling. As the coupling increases, this resonance 
splits forming a doublet for which the quantum numbers will be J n = | + and | + . The first 
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TABLE III: Quantum numbers of allowed even-parity basis states (J 7 * < | ). 
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of these we identify with the state N = 5 (responsible for the resonant behavior), coupled 
with an elastic background arising from the effect of the N = 4 entry. The second of the 
doublet we link to the state N = 10 and expect that it is attached to a background state 
for which N = 9. As the deformation (l/^l) increases, the splitting of the doublet, and the 
resonance widths, increase. Also all of the J-components contribute since entries N=4, 5, 6, 
7, and 8 are involved for J = | and those of N=9, 10, 11, 12, 13, and 14 are so for J = |. 

Hence the bound state \ + is expected to generate a doublet of compound resonances, and 

by the same mechanism, the bound state | + when coupled to the 2 + state should generate 

a quintuplet with spin-parities J = | + , | + , | + , | + , and | + . 

To verify this scheme, we have made calculations starting with the physical situation of 
having /3 2 = —0.52, and then form cross sections from calculations in which the coupling is 
gradually reduced. By that means it is possible to track each state/resonance continuously 
and so identify the underlying base nature of each state. It is assumed that each resonance 
conserves its identity in the adiabatic transition from the physical to the unperturbed limit. 
We portray the results in Fig HJ which shows the general trend that the background cross- 
section values near threshold increase in size as l/^l decreases. That is consistent with the 
strong sub-threshold s-wave | state moving closer to threshold as \/3z\ decreases. Then 

what was just a sub-threshold state, the | in the physical case calculation, moves into 
the positive energy regime with decrease in |/? 2 |. That state has an extremely small width 
(< 10~ n MeV) so that while its existence is known, its strength has not been ascertained. 
The actual resonances also have unique trends with decreasing \j3 2 \. With their identification 
(the entry 'n' specified in Table ITT|) . the narrow | (n = 10) and the first | (n = 6) track 
smoothly with decreasing width till they have essentially the same centroid energy and 
vanishing width; features that we discuss subsequently. The very narrow | + (n = 12) 

resonance initially sits upon the | + resonance shape and its centroid increases in value with 
decreasing |/3 2 |. It remains a feature upon the | + shape until l/^l ri 0.35, after which that 
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FIG. 4: (Color online) Total elastic cross sections for n+ 12 C scattering as functions of neutron 
energy showing the effects of reducing the value of j3 2 . Details are as described in the text. 

| + centroid moves up in energy to be one of the quintuplet of sharp resonances at ~ 4.3 

MeV. The other four members of the quintuplet are the very broad second | + (n = 7), the 

reasonably broad strong | + (n = 11), the sharp | (n = 13), and the reasonably broad but 

weak | + (~4.6 MeV) resonances. The last is not readily seen until, with decreasing |/3 2 |, the 

width of the nearby | + resonance decreases sufficiently. In the figure, this | + resonance is 
seen as the highest energy track that curves back in energy to meet the rest of the quintuplet 
when (3 2 — 0- The | + resonance sits initially on the tail of the | + resonance appearing to 

gradually decrease in size ( because the | + resonance is contracting as \0 2 \ decreases) before 
some enhancement at small values of \ft 2 \ due, then, to the proximity of all members of the 
quintuplet. The centroids of the two broad resonances (the | + and the | + ) of this quintuplet 
cross as \/3 2 \ decreases and while each width gradually decreases, they combine to form quite 
distinctively changing cross sections until, as (3 2 — > 0, they become very narrow and tend to 
the same centroid energy. Of course when f3 2 is exactly zero, since the widths are vanishingly 
small, a calculated cross section loses almost all trace of individual compound resonances. 
Nevertheless by means of the RI procedure, the resonance centroid energies may still be 
(and were) found. 
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FIG. 5: (Color online) Variation of energies of the \ bound and doublet, and of the | bound 
and quintuplet of states with deformation. The numbers in brackets attached to each curve are 
the specific values of 2J. 

In Fig. |5] we show our results of the calculations of the centroids of the positive parity 
resonances and sub-threshold bound states. Consider first the states built upon a | + single 

particle state. They are the quintuplet of resonances having spin-parities J n = | + , | + , | + , 

| + , and | + plus the sub-threshold (bound) state in 13 C. Their energy values when (3 2 = 
(i.e. ~ —0.15 and 4.28) equate to the target excitation 4.28 — (—0.15) = 4.43 ~ e 2 , and 
so we conclude that the quintuplet indeed is a set of compound resonances due to coupling 
of a | + neutron with the 2 + excited state of 12 C and that the Ods nucleon binding to the 

ground state of 12 C is ~ —0.15 MeV. In a similar way, the | bound state (n = 3) has a 
binding to the ground state of 12 C of ~ —0.83 MeV and the doublet of resonances (n = 6 
and 10) relate in energy to identify their origins as the | , | doublet of states generated 

by that | + neutron being coupled to the 2 + excited state in 12 C. In the zero deformation 
limit, the energies of the bound and collapsed doublet states are ~ —0.83 and 3.60 MeV 
thereby differing by the 4.43 ~ e 2 value. 

This leads to a complete understanding of the spectroscopy of the even-parity states in 
13 C as they are revealed by low energy n+ 12 C scattering. Note that the two | + states (n = 6 
and 7) that come from analysis solely of the scattering data, basically are indistinguishable 
because only J n is conserved. Nevertheless we expect that n = 6 is part of the doublet and 
n = 7 is part of the quintuplet because of the continuity in the limit (3 2 — > 0. Of course, 
similar information should be obtained by large-space structure model analyses of the wave 
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TABLE IV: Properties of the even-parity quintuplet of resonances in the limit (3% = 0. 
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functions in the two cases. Work is in progress to do just that. 

We note that the = limit values of energies of the even-parity states discussed tend 
only approximately to the same values. There are small but significant shifts. That effect 
is due to the residual operator character of the interaction potentials. While the coupling 
between the channels due to deformation was removed, there remains some linkages arising 
from the diagonal potentials being 

V e (r) = V cc (r) = [Vo + £{£ + l)Vu) fo(r) - ^W ls A{r) [1 • s] cc + V„f (r) [s ■ I] cc . (2) 

This is the reduced form of Eq. (B.5) of Ref [H, in the limit @2 — > 0, where it is meant 
that the channel index c denotes the set of quantum numbers {J 71 " , I,j,£}. Clearly the 
above discussion fails to hold exactly in presence of a spin-spin interaction. How much so 
is illustrated by the entries in Table IIVI In this table we have selected the states of the 
even-parity quintuplet, and show in column 3, the calculated values of the unperturbed 
energies. They reveal the small defects of convergence caused by the spin-spin (target) 
specific potential term. It is easy to recognize that the members of the even-parity quintuplet 
(in the unperturbed conditions) are pure states of the type N = 3, 7, 12, 17, and 22, as given 
in Table ITTT1 with quantum numbers c = 2, §, 2}, {| + , 2, §, 2}, {| + , 2, §, 2}, {| + , 2, §, 2}, 
and {| + ,2,|,2} respectively. These states differ (only) in the total angular momentum J 
and therefore are discriminated (only) by the spin-spin part of the potential. 

In particular, from Eqs. (B8), (B9) and (B10) of Ref. pj, it may be seen that the spin- 
spin potential eigenvalues depend on the quantum number J through functions K(£,j,j', J). 
Those numbers are listed in column 4 of Table IIVI To verify the proportionality of the 
values of E th and K, we give the differences AE th between adjacent states in the list (i.e. 
£?th( n — 7) — -Eth(n = 4), etc.) in column 5. Then, in column 6, the associated differences 
AK (namely K(n = 7) — K(n = 4), etc.) are given while in column 7 we present the ratios 
■xw~- The value reported in the last column is constant within 0.4 percent, most readily 
confirming the assumption that the degeneracy of the multiplet in the unperturbed limit 
is broken only by the spin-spin interaction. Therefore we have repeated calculations with 
$2 = but additionally with the spin-spin potential set to zero. All other parameters remain 
with those values given in Table HJ The energies of the states then found are presented in 
Table IV1 Therein also, in the rightmost column, we give the gap energies AE n which are the 
differences between the resonance energies of the listed value and of the bound state energy 
of the particle coupled to the 12 C ground state. Therein the states for n = 1 , 3 , and 9 are 
taken as coupling a Opi, lsi, and a Ods neutron to the ground state of 12 C. Then all of the 
remaining differences shown in the last column (except for the n = 2 case) are exactly the 
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TABLE V: 


Sub-threshold bound states and resonances in n+ 12 C when {5% = and V ss 


= 0. 


n 






AE n 


1 


l- 

2 


-4.7017 





2 


1- 

2 


2.9525 


7.6542 


3 


1 + 

2 


-0.8376 





4 


1 + 

2 


4.2839 


4.4389 


5 


3 _ 
2 


-0.2627 


4.4390 


6 


3 + 
2 


3.6013 


4.4389 


7 


3 + 
2 


4.2839 


4.4389 


8 


5~ 
2 


-0.2627 


4.4390 


9 


5 + 
2 


-0.1550 





10 


5 + 
2 


3.6013 


4.4389 


11 


5 + 
2 


4.2839 


4.4389 


12 


7 + 
2 


4.2839 


4.4389 


13 


9 + 
2 


4.2839 


4.4389 



assumed excitation energy of the 2 + state in 12 C (e 2 ). As already observed, our assumption is 
that the even-parity quintuplet, n = 4, 7, 11, 12, and 13, derives from the bound state giving 
n = 9 (and therefore the definitions of the gap energies are AE^ — E4 — Eg } /\Ej — Ej — Eg 
etc), while the even-parity doublet, n = 6 and 10, derives from the bound state that lead to 
n = 3 (and therefore AE 6 = E e — E 3 , AE W = E 10 — E 3 ). The n = 2 entry results from the 
coupling of the Opi (that gave n = 1) to the third state we chose to include in the target 

spectrum used in MCAS. This is the only effect of the excitation of the excited state in 
12 C at least for the energy range we have considered. To explore more effects of the 0^" state, 
higher energy results need to be analyzed. But then one would also need consider effects of 
other states such as the strong collective 3~ (9.63 MeV). 

As far as the other odd-parity states are concerned, the values listed in Table |V| strongly 
suggest that the states n = 5 and n = 8 (which is reproduced as a near-threshold bound 
state in our calculations) are the result of coupling a Opi neutron (as identified in the n = 1 

level) with the 2 + excited state in 12 C. All these conclusions about the odd states may be 
drawn only after dropping the spin-spin term, because the spin-spin interaction is stronger 
for the odd states (see Table HJ) and masks completely the effect of the convergence in the 
unperturbed limit. 

Finally with the n+ 12 C system, in Fig. |U] we show the even-parity components of the 
cross sections but now when a small finite deformation = —0.2 has been used. In this 
case only cross sections in the energy range 2.4 to 5.5 MeV are displayed which suffice to 
show three groupings of interest. The resonances seen in the total elastic cross section for 
this (small) deformation in sequence as energy increases have spin-parities, | , | + , | + , | + , 

| + , | + , | + , and | + . The total cross section is depicted by the solid curve. The notation for 
the separate component even-parity cross sections are as given for Fig. El By comparison 
with results shown previously, we note that: 

1. The resonances can be assembled into 3 groups: an odd-parity singlet, an even-parity 
doublet and an even-parity quintuplet. 
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FIG. 6: (Color online) Even parity component contributions to the elastic n+ 12 C cross sections 

I + and |" 



in the case of small deformation, f5 2 = —0.2. In this case the narrow X and § resonances lie at 



(lab) energies of 4.0 and 4.8 MeV. 



2. The | resonance (n = 4), though registered through the RI procedure but which 
was not evident as a peak in the cross section in Fig. Q now is seen as a little bump 
in the cross section above 5 MeV. Of the even-parity components, the | + one which 
dominates the lowest energy values, clearly identifies that highest energy resonance in 
this figure as having that spin-parity. 

3. The widths are smaller than those found in Fig. ^ and as was also very evident in 
the 3-D diagram (Fig. HJ). They all tend to zero as \(3 2 \ is gradually taken to zero, so 
confirming the nature of all of the compound resonances analyzed. To illustrate that 
variation in the widths T, a select set of those for some of the larger resonances are 
plotted versus the deformation f3 2 in Fig.0 As foreseen by theory the widths tend 
to zero as |/3 2 | 2 . Note that in one case (n = 7), T th (J n = | + ) initially increases as \fl 2 \ 
decreases. This is due to the coupling being quite strong so that the physical regime 
is far from the unperturbed limit. 

4. We have verified that with the small value for deformation, the n = 8 | state is a 
resonance, a bound state in the continuum. For the physical value of (3 2 = —0.52, it 
is an actual sub-threshold bound state. 
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FIG. 7: (Color online) Deformation dependence of widths, r t j 1 (J ?r ) (in keV), for some of the 
(broader) compound resonances. Those from cases with n = 4, 6, 7, 10, and 11 are portrayed by 
the long dashed, dot-dashed, solid, thin solid, and dashed curves respectively. 

III. THE 13 N SYSTEM 

In analyzing the 13 N system we have assumed charge symmetry, and so use the parameter 
values given in Table H] for the n+ 12 C system to evaluate p+ 12 C scattering and to assess 
properties of the sub-threshold compound nucleus 13 N. Only a Coulomb interaction has 
been added. For simplicity we have chosen the Coulomb potential to be that of a uniformly 
charged sphere so that we have only one new parameter, the Coulomb radius R c . The best 
agreement with the experimental data has been obtained with the value R c = 2.4 fm. 

The results, centroid energies and widths, are listed in Table |VT] and compared therein 
with experimental data |8j. The comparison is quite good with the mean square error /i, as 
specified by Eq. (jTj), relative to 13 states of Table IVT1 being 0.2563 MeV. The entry in row 
11 has been excluded from the sum since there is no experimental counterpart to that state. 
The average of the differences between pairs of states in 13 N and 13 C with reference to the 
experimental and theoretical spectra respectively, are the following: 

(£ exp ( 13 N)-£ cxp ( 13 C)) = 3.01 MeV ; (i?th( 13 N) — i? t h( 13 C)) = 2.54 MeV . (3) 

Thus, through coupling with the 2 + level, the bound state underlying the n = 1 entry gen- 
erates the compound odd-parity doublet n = 5 and 8 in the tabulation, and that underlying 
the resonance for n = 3 generates the quasi-compound even-parity doublet n = 6 and 10. 
Likewise, the resonance n = 9 links to the quasi-compound even-parity quintuplet formed 
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TABLE VI: Comparison between experiment and theory, p+ C. 
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5 + 


1.6035 


23.5 


0.6379 


0.899 


10 


1+ 


4.4205 


5.5 


4.1794 


8.86 


11 


1+ 






6.5281 


726 


12 


?+ 


5.2115 


4.5 


5.1234 


1.50 


13 


1+ 
2 


7.0565 


140 


6.8341 


153 


14 


5 + 
2 


9.5865 


215 


9.7895 


910 



TABLE VII: Sub-threshold bound states and resonances in p+ C when 02 = and V ss = 0. 



n 




£th (MeV) 


AE n 


|r th (k. 


1 


l- 

2 


-1.6785 







2 


1- 
2 


5.9757 


7.6542 





3 


1 + 
2 


1.0217 





468 


4 


1 + 
2 


6.6841 


4.4389 


115 


5 


3~ 
2 


2.7604 


4.4389 





6 


3 + 
2 


5.4606 


4.4389 


467 


7 


3 + 
2 


6.6841 


4.4389 


115 


8 


5~ 
2 


2.7603 


4.4388 





9 


5 + 
2 


2.2452 





115 


10 


5 + 
2 


5.4606 


4.4389 


466 


11 


5 + 
2 


6.6841 


4.4389 


115 


12 


7 + 
2 


6.6841 


4.4389 


115 


13 


9 + 
2 


6.6841 


4.4389 


115 


14 


5 + 
2 


9.8995 


7.6543 


115 



by the n = 4, 7, 11, 12, and 13 entries in Table IVll Finally, the n = 1 bound state when 
coupled with the second + level gives rise to the single compound resonance for n = 2. 

For confirmation, in Table IvTTl we give the results of calculations of the p+ 12 C system 
in the unperturbed limit (/?2 = 0) with the spin-spin term set to zero. The fourth column 
shows the differences which are consistent with the couplings we suggest. All the members 
of each multiplet tend to the same limit, and the difference between this limit and the 
relative generator single particle state corresponds to the energy of the relevant target excited 
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state involved in the coupling. But there are some differences. For example, the proton 
calculations predict four | + resonances in the energy range selected. That is one more 

than found from the neutron calculations in an equivalent energy range. Also, while the | + 
resonances n = 9, 10 and 11 are the single particle state, a component of the even-parity 
doublet, and a component of the even-parity quintuplet as found in the neutron case, the 
n = 14 resonance is new. Our fa — > analysis indicates that this (extra) state originates 
from the coupling to the excited 0^~ state in 12 C. 

An analysis of the widths listed in Table IVHI is interesting. In the past |9( , and derived 
from a purely mathematical example, it was noted that the width of the compound resonance 
approaches zero as the coupling approaches zero (because its origin was a bound state), while 
the width of the quasi-compound ones tend to the natural width of the resonance from which 
such originated. It is evident in Table rvTTl that the odd-parity doublet (n = 5 and 8) and 
the odd-parity resonance (n = 2), are compound resonances as their widths vanish in the 
limit of zero deformation. The even-parity quintuplet (n = 4, 7, 11, 12, and 13) and the 
resonance (n = 14) have the width of the common resonance partner (n = 9) and so are 
quasi- compound resonances. The width of the quasi-compound even-parity doublet (n = 6 
and 10) is that of the partner resonance for which n = 3. 

Thus there is a one to one correspondence between the states calculated by the model 
and those experimentally known, with the exception of the | + resonance (n = 14). But that 
resonance, while missing in n+ 12 C experimental data, has an observed partner in p+ 12 C 
system. Also, while it is clear from the comparisons to be made within Table IVII that the 
calculated energy centroids generally are in good agreement with data and that the widths 
concur in order of magnitude, there are some exceptions. In particular, the | (n = 8) state 
is not well reproduced. Also state n = 3 is very close to threshold and, in spite of the small 
difference between theory and experiment, calculations make it be a sub-threshold bound 
state instead of a resonance. 

Of course we have adopted a strict charge symmetry assumption and have used a very 
simple form for the Coulomb interaction. It should be possible to find improved 13 N prop- 
erties from the MCAS approach with but slight changes of these potentials. 

A. Proton scattering from 12 C; cross section and analyzing power 

Since Coulomb scattering amplitudes diverge at zero degree scattering, measurements of 
proton scattering do not lead to total elastic scattering cross sections. Instead the usual 
procedure is to find cross sections at fixed scattering angles and/or differential cross sections 
at fixed energies. However, as the MCAS approach yields complete scattering (S-)matrices, 
such angular observables are readily predicted. 

For energies to 7 MeV, proton elastic scattering from 12 C at fixed (cm) scattering angles 
of 54° and at 90° are shown in the top and bottom panels of Fig. |H] respectively. Therein 
our calculated cross sections, shown by the continuous lines, are compared with available 
data. With energies in the laboratory frame, the data measured at 54° were taken from 
Ref. ^(j (open circles) while those depicted by filled black circles come from Ref. Jn|]. At 
54°, the |" (n = 5) resonance in the data is very well reproduced. The other resonances, 

| + (n = 10), | + (n = 6), and the broad | + (n = 7), also are well reproduced in shape and 
require but a small shift in energy. That need for a small energy shift is also evident in their 
centroid values that are listed in Table IVII The sharp peak shown in the theoretical curves 
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FIG. 8: (Color online) Cross sections from proton elastic scattering from 12 C for two (cm) scat- 
tering angles, 54° (top panel), 90° (bottom panel). Details are given in the text. 



near 3.2 MeV is the | (n=8) resonance. This state has been experimentally detected at 
higher energies, These findings are confirmed by comparison of the calculated results with 
the data taken at 90° and displayed in the bottom panel. 

Usually analyzing powers A y from nucleon-nucleus scattering are more sensitive to details 
of structure than are the differential cross sections. Two results for this observable are 
presented in Fig. In the top panel of this figure, analyzing powers for proton elastic 
scattering at 90° (cm) and for (lab) energies 1 to 8 MeV are depicted. In the bottom panel, 
angular variation of the analyzing power for a fixed energy of 3.5 MeV is shown. The 
experimental data lj| are depicted by open circles and our calculated results are given by 
the solid curve. The resonances seen in this figure are those that were discussed in relation 
to Fig. |HJ The bottom panel gives an angular variation of analyzing powers for an energy 
(3.5 MeV) that is well removed from any strong resonance influence. The calculated result 
then reflects what the model predicts as background effects. The agreement with data is 
good. Note that no additional parameter variation was done to achieve this agreement. 

To investigate deformation effects on differential cross sections, we calculated them again 
using the values of parameters as listed in Table but with the deformation parameter 
$2 = —0.2. The results are shown in Fig. El for a (cm) scattering angle of 90° and (lab) 
energies ranging between 1 and 6 MeV. This figure reveals how the single particle and com- 



pound resonances behave with decreased deformation. The | (n = 3) and 



5 + 



9) 



single particle resonances (the first of which was not seen in the preceding figure because it 
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FIG. 9: (Color online) Experimental analyzing powers from proton scattering from 12 C compared 
with our theoretical calculations. Data and results taken for a range of energies at a fixed scattering 
angle of 90° are shown in the top panel, while angular variations at a fixed energy (of 3.5 MeV) 
are shown in the bottom panel. 



lies too low in energy) maintain their natural width, while the | (n = 5) and | (n = 8) 
compound resonances sharpen. With fa — —0.2, already they are very narrow. The small 
bump between 5 and 6 MeV corresponds to the | + (n = 10) resonance which maintains its 
natural width but, being quasi-compound, is strongly damped as the unperturbed limit is ap- 
proached. Thus the behavior of the three kinds of resonances (single particle, compound and 
quasi-compound), as the coupling constant is decreased is just that foreseen previously 0. 
Namely, as f3 approaches zero, single particle resonances conserve their shapes, compound 
resonances reduce their widths gradually keeping constant heights, while quasi-compound 
resonances reduce their heights gradually keeping constant widths. 



IV. CONCLUSIONS 

We have analyzed the low energy spectra of 13 C and 13 N and low energy data from the 
elastic scattering of neutrons and of protons from 12 C. Our method of analysis was that of a 
multi-channel algebraic scattering theory with which both the (positive energy) scattering 
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FIG. 10: Energy variation of the theoretical differential cross section for proton scattering from 
12 C at a cm scattering angle of 90°. The calculation was made using = —0.2. 



of nucleons from 12 C and (by using negative energies) the sub-threshold bound states of 
the compound nuclei could be predicted. The method, based upon sturmian expansions of 
a matrix of interaction potentials, ensures that all sub-threshold as well as resonant states 
within the chosen range of energies are found. Also the approach can be adapted so that the 
Pauli principle is not violated even when a collective model is used to define those interaction 
potentials. In the cases studied, just such a rotational model was used for that purpose and 
the 0+ (ground), 2+ (4.438 MeV), and the Oj (7.96 MeV) states in 12 C taken as active. 

The results of our analysis are summarized in two figures, Fig. ^2 and for the com- 
pound nuclei 13 C and 13 N respectively. Each displays three ladder diagrams; the first two are 
the results of our calculations (unperturbed resulting when (3 2 = 0) and they are compared 
with the experimental values (far right). Considering Fig. for 13 C first, we note that it 
supports three (dominantly single particle) bound states, as is clear from the unperturbed 
(@2 = 0) spectrum depicted on the left of the figure. Those states are highlighted in that 
spectrum as dark lines. In addition the coupling of the incoming neutron with the excited 
levels in the target 12 C then gives rise to meta-stable states whose presumed unperturbed 
{02 — 0) configurations are represented by lighter lines. They are connected with what we 
deduced as a partner by the curved lines. The energy gaps between these are e 2 and e 3 
as relevant for the two possible cases. These unperturbed resonances are infinitely narrow 
(r t h — > 0) and fully degenerate when one neglects the small effect of the spin-spin interac- 
tion in this representation. The finite deformation (/% = —0.52) splits these components to 
yield the predicted resonances as shown in the middle of the set and which compares very 
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FIG. 11: (Color online) The genesis of the theoretical spectrum for 13 C schematically represented 
and compared with the experimental spectrum. 



favorably with the experimental spectrum given in the box to the right. Similar conclusions 
can be drawn from the spectra for 13 N that are shown in Fig. El Because of the Coulomb 
energy shift, there is only one sub-threshold bound state, and the MCAS approach assuming 
charge symmetry gives just that and with the correct spin-parity. As with the 13 C system, 
accounting for the Pauli principle is crucial. Otherwise numerous spurious levels result. An 
interesting difference however is that the Coulomb effect transforms the | + and | + bound 
states to be single particle resonances, and their products due to coupling with excited core 
states become quasi-compound resonances. The different behavior of compound and quasi- 
compound resonances in the limit 02 — > confirms statements formulated in a previous 
paper 

We conclude noting that, with a unique set of potential parameters the MCAS approach 
reproduces data. The results facilitate interpretation of the phenomenology of both the 13 C 
and 13 N systems in the considered energy range in both a satisfactory and a fairly exhaustive 
way. 
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